A MULTIPARAMETER SUMMATION FORMULA 
FOR RIEMANN THETA FUNCTIONS 
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Abstract. We generalize Warnaar's elliptic extension of a Macdonald multi- 
parameter summation formula to Riemann surfaces of arbitrary genus. 



We start by a brief outline of the general theory of hypergeometric type series 
built from Jacobi theta functions ^3 (its extension to integrals ^Sj is not 
touched at all). Within this approach, univariate elliptic hypergeometric series are 
defined as formal series Y] c n for which h(n) = c n+ \/c n is an elliptic function of 
n considered as a continuous complex variable. Normalizing cn = 1, we see that all 
coefficients c„ are obtained as products of h(k) or l/h(k) for different k G Z. 

Any elliptic function of order r + 1 can be represented as 

_ ^(up + n, . . . ,u r + n; <t,t) = ^ 9(t Q q n , ■ ■ .,t r q n ;p) 
0i(t> o +n, ...,v r + n;cr,T) 9(w q n , . . . ,w r q n -p)' 

k k 

9i(u , . . . , u k ; a, t) = Jj0i(Ki;o-,r), 9(t , . . -,t k ;p) = J|6»(^;p), 

1=0 2 = 

with a free variable z G C, the base p = e 2vlT satisfying the constraint \p\ < 1 (or 
Im(r) > 0), and arbitrary q = e 27ria . The function 

0i(u;«7,t) = -i jr (-l)V a +D 2 /8 ? (W/2)« 

k— — oo 

= ip 1 ^q- u/2 ( P ;p) oo 0(q u ;p), u G C, 
is the standard Jacobi #i-function. The short theta function 9(a;p) has the form 

oo 

Q(a;p) = (a;p)oo{pa^ 1 ;p)oo, (a;p)oo = JJ (1 - ap n ), (2) 

71 = 

and obeys the properties 

9(pa;p) = 9(a- 1 ;p) = -a- 1 e(a;p). (3) 

We use also the convenient elliptic numbers notation: 

[it; cr, r] = 6\(u; cr, r), or [u] = 6i(u). 

The parameters Ui and Vi in are defined modulo shifts by cr -1 , which will 
not be indicated in the formulas below. They are connected to ti and Wi as 

U = q Ul , Wi = 
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and satisfy the constraint 

r r r 

^2{ Ui -vi)=o, or n* i= n^' ( 4 ) 

i=0 i=0 i=0 

which guarantees double periodicity of the meromorphic function h(n): 
h(n + o" -1 ) = h(n), h[n + rcr -1 ) = h{n). 
For the entire function [u] we have [— u] = —[u], 



\u + a- 1 ] = -\u], \u + ra- 1 ] = - e -™-2™u M 



(5) 



and 



[u;a,T + l}=e^ i [u;a,r}, (6) 
[u;a/r,-l/r] = -i{-i T f/ 2 e™ 2 « 2 /r [u . ^ t]j (7) 

where the sign of (— ir) 1 / 2 is fixed from the positivity of its real part condition. 
The latter two transformations generate the modular PSL(2, Z)-group, 

ar + b a 

T ^ — Tl> ° — 71' ( 8 ) 
ct + a ct + a 

where a,b, c, d £E Z and ad — be = 1. 

For the unilateral series X)neN c ™' we conventionally normalize wo — q (or 

t>o = 1). Then, the elliptic function h(n) generates the single variable elliptic 

hypergeometric series: 

„ ( to,...,U \ \^ 9(to,ti, . . . ,t r ;p; q) n 

r+iE r [ '>Q>P'' z ) = l^Tr r z . (9) 

\wi,...,w r J ^ 6{q,wi,...,w r ;p;q) n 

The elliptic shifted factorials are defined as 

k n — 1 

6(t ,...,t k ;p;q) n = Y[ Y\_0(tmq J ;p), 

m=0 j=0 

or in the additive form: 

k n—1 

[u Q , Uk]n \ \[ [u m + J*! °i A- 

m=0 j=0 

A further generalization of these functions is constructed from h(n) equal to an 
arbitrary meromorphic Jacobi form in the sense of Eichler and Zagier 5 . The 
corresponding series were called in |f3J as theta hypergeometric series. In this 
classification, theta hypergeometric series (0 are called balanced if their parameters 
satisfy the constraint (@J. 

Elliptic hypergeometric series © are called totally elliptic if h(n) is also an 
elliptic function of all free parameters entering its #i-functions arguments. It is not 
possible to have this property if all Uj and Vi are independent. Therefore such a 
requirement results in the following additional constraints |13j : 

uq + 1 = ui + v i = . . . = u r + v r , or qto = t\Wi = . . . = t r w r . (10) 

These relations are known as the well-poisedness condition for plain and basic hy- 
pergeometric series Thus the total ellipticity concept sheds some light on the 
origin of the restrictions l|l(J|l . Totally elliptic hypergeometric series are automati- 
cally modular invariant. 
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The balancing condition for well-poised series contains a sign ambiguity: ti ■ ■ ■ t r 
= ±g( r+1 )/ 2 £g r ~ 1 ' ) / 2 . It clearly shows also distinguished character of the parameter 
to- It is convenient to express t r in terms of other parameters. In this case the 
function h(n) is invariant with respect to the independent p-shifts to — > p 2 to and 
tj — > ptj, j = 1, . . . ,r — 1. However, for odd r = 2m + 1 we can reach the to - * pto 
invariance provided we resolve the ambiguity in the balancing condition in favor 
of the form t\ ■ ■ ■ t^m+x = +<z m+1 t n , which is precisely the condition necessary for 
obtaining nontrivial hypergeometric identities. Thus the notion of total elliptic- 
ity (with equal periods) uniquely distinguishes the correct form of the balancing 
condition. 

The elliptic analog of the very-well-poisedness condition consists in adding to 
(|10|1 of four constraints 

. ,1/2 . ,1/2 1/2 _i/2 . ,1/2 1/2 ri1 s 

t,._ 3 = t ' q, t r _ 2 = -t <7, U-i = V qp ' , t r = -t ' qp 1 . (11) 

After simplifications, the very-well-poised theta hypergeometric series r +\E r take 
the form (no balancing condition is assumed) 



-4 



In the limit p — > 0, functions (|12|) are reduced to the very- well-poised basic hy- 
pergeometric series r _i<^ r _2 of the argument — qz 8 . For even r the balancing 
condition for O has the form t x ■ ■ ■ t r _ 4 = ±q ( ~ r - 7 ^ 2 t^'~ 5)/2 . For odd r = 2m + 1, 
it is appropriate to call (|12f> balanced if 

2m -3 

n h = r-x- 2 - 

Then the function h(n) is elliptic in uo, i.e. it is invariant with respect to the shift 
to — > pto- Conditions l|ll|) do not spoil this property because the change to — > pto 
is equivalent (due to the permutational invariance) to a replacement of parameters 
(|llfl by t r _3,pt,._2,pir-i, tr an d, so, it does not have an effect upon h(n). The 
modular invariance (which is not immediately evident due to the dependence of 
parameters in (|11|) on the base p) is preserved due to a similar reasoning. 

For the first time, the very-well-poised elliptic hypergeometric series (with a 
different way of counting the number of parameters and the choice z = — 1) have 
been considered by Frenkel and Turaev [7]- Their work was inspired by exactly 
solvable statistical mechanics models built by Date et al In an independent 
setting, such functions have been derived by Zhedanov and the author by solving 
a three term recurrence relation for a self-similar family of biorthogonal rational 
functions [T7|. 

As shown by Frenkel and Turaev [7], the following summation formula is true: 
■A 6(t Q q 2k ;p) -A- 9{t m ;p; q) k k 

l<r<s<3 @{Qbo/trts',P) <?)n . , 

9{qt /tit 2 h;p- 1 q) n n;Li 8(qto/U;p; q) n ' 

where Ili=i U = qt 2 , and t 4 = q~ n , n E N. This means that the special terminating 
very- well-poised balanced lo^g series with z — — 1 is summable. If we take the 



4 



VYACHESLAV P. SPIRIDONOV 



limit £3 — > to/t\ti, then the termination condition cancels out and (|13J) becomes a 
summation formula for an indeterminate %Ej series 

0(toq 2k ;p) A 0(t m ;p;q)k ^ fe _ 0{qt ,qt 1 ,qt 2 , qh/hhip; q) n 
k~o ^* o;p ) m=b s (i t o/ t m;p;q)k d(q,qto/t 1 ,qt /t 2 ,qtit 2 ;p- 1 q)n' 

For p — > 0, we obtain an indeterminate summation formula for a very-well-poised 
balanced 6^5 basic hypergeometric series. In the additive notation, equality l|14|) 
takes the form 

\[u Q + 2k] [u ,u 1 ,u 2 ,u -u 1 -U2]k 



E 

fc=0 



[Uo] [1, 1 + 7X0 - Til, 1+ Uo - 7X2, 1 + 7X1 + U 2 ] k 

_ [1 + TXq, 1 + Ml, 1 + 7X 2 , 1 + TXp - Ui - U 2 ] n 
[1, 1 + 7X0 - 7X1, 1 + 7X0 — 7X2, 1 + 7X1 + U 2 ] n 



(15) 



The general theory of multiple elliptic hypergeometric series is much less de- 
veloped. Several examples of such series associated with root systems have been 
investigated by Warnaar [T5]. van Diejen and the author 0J^||E3, Rosengren [T2] . 
Kajihara and Noumi [S], Rains with the latest contributions coming from the 
Gustafson's talk at the present conference. 

We expect that the approach of [21 to elliptic hypergeometric series is generaliz- 
able to multidimensional Ricmann theta functions associated with algebraic curves. 
In this case, we can take as h(n) a meromorphic function upon a compact Riemann 
surface S of an arbitrary genus g. However, the things are much less transparent in 
this situation. In [5], we started to discuss such generalizations and suggested an 
extension of sum (|15|) to a series built from the Riemann theta functions. However, 
there was a flaw in the formulation of this result. The main purpose of the present 
note is to lift the Macdonald multiparameter sum described in Theorem 2.27 of 1 
to Riemann surfaces of arbitrary genus and to correct a related statement of |15j . 

The Riemann theta function of g variables TXi, . . . , u g € C with characteristics 
a = (ai, . . . , a g ) and j3 = . . . , (3 g ) is defined by the g-fold series of the form 

9 

6 a:/3 (u; fi) = ^2 exp{7ri ^ ( n j + ^j)^jk(rik + £*fc) 

nGZs j,k=l 

9 

+ 2m^(7x; + f3 3 ){n 3 + a,)}, (16) 
i=i 

where Qjk is a symmetric matrix of periods. For general abelian varieties the matrix 
fljk is arbitrary, and for the Jacobians associated with the Riemann surfaces, it is 
generated by a basis of holomorphic differentials lj = (wi, . . . ,u g ) (see, e.g., [TU)'). 
It is evident that for arbitrary j we have 

6^(7X1, . . . , Uj + 1, . . . ,tx s ; n) = e 2 ^e Q ,0(u; Q). (17) 

Analogously, we have 

6^(7X1 + fill,, ••••», ■ fi 9 k;fi) e -^"-2-i(A+«* t )e a)/3 ( U ; Cl), (18) 



where k = 1, . . . , g. 
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We denote as IY2 a subgroup of the symplectic modular group Sp(2g, Z) gener- 
ated by the matrices 

l=( a c h d ) eSp(2g,Z), 

such that diag(a t b) — diag(c*d) = 0mod2. The action of this group upon the 
matrix of periods Q and the arguments u of the theta function is defined as 

n' = (an + ^icn + d)- 1 , u' = u*(on + &) -1 . (19) 

Analogously, we define the tranformed characteristics 

( a' \ _( d -c\fa\ If diag^d) \ 
\ J ~ \ -b a J\P J + 2 \ diag^d) ) ' 

Then, the 1Y2 group transformation law for theta functions has the form 

G a ,, ,{u'; n') = Cv / dit(dY+d)e' riut ( cn+d )" lcu e a ^(u; fi), (20) 

where £ is an eighth root of unity 
We denote as 

,6 

Vj(a,b) — / uij, a,b E S, (21) 

J a 

abelian integrals of the first kind. The characteristics ctj,f3j G (0, 1/2), such that 
4:^2 j— 1 a jPj — 1 (mod 2), are called odd. We denote theta functions with arbitrary 
(nonsingular) odd characteristics as [u], u S C 9 , and take the convention that 
[ui, . . . , Ufc] = nj=i[ u j]- for such functions, we have [— u] = — [u]. It is known that 
theta functions on Ricmann surfaces satisfy the Fay identity 0: 

[u + v(a, c), u + v(6, d), v(c, b), v(a, d)} 

+ [u + v(6, c), u + v(a, d), v(a, c), v(6, d)] (22) 
= [u, u + v(a, c) + v(6, d), v(c, d), v(a, 6)], 

valid for arbitrary ugC 9 and a,b,c,d G 5* (note that v(a, c) + v(6, d) = v(6, c) + 
v(a, d)). Here we can replace theta functions independent on the variable u by 
prime forms E(x, y) since the cross ratios of these theta functions coincides with 
the cross ratios of the corresponding prime forms (and therefore they are indepen- 
dent on the choice of odd characteristic for theta functions). However, we find it 
more convenient to work with one theta function [u]. If we consider the function 
[u + v(a, 6)], then, in the appropriate normalization, transformations (I17|) and 118|) 
correspond to cyclic moves of the point b (or a) on the Riemann surface along the 
Ai and Bi contours, such that f. u)j — #y and j B ujj = Clij, respectively |l()j . 

We would like to describe a summation formula for a particular series built from 
the Riemann theta functions. 

Theorem 1. We take a nonnegative integer n and consider n+1 arbitrary variables 
Zfc e C 9 and 4n + 4 different points on a Riemann surface dfe, 6fe, Cfe , dk G S, k = 
0, . . . , n. Then the following multiparameter summation formula for Riemann theta 
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functions on algebraic curves takes place 

n 

^2[z k + v(b k ,c k ),z k + \r(a k ,d k ),v(a k ,c k ),v(b k ,d k )} 
k=o 

k-l 

x Y[[zj,Zj +v(a j ,c j )+v{b j ,d j ),v(c j ,d j ),v(a j ,b j )] 

71 

x II t z J + v ( a i> c i)' z J + v ( & j = d j)> v ( c j> & j)> v ( a J> d i)] ( 23 ) 
i=k+i 

n 

= JJ[z fe ,z fc + v(a fc ,c fc ) + v{b k ,d k ),v(c k ,d k ),v(a k ,b k )} 

fe=0 

n 

- J]jz fe + v(a fe ,c fc ),z fc + v(& fe ,d fe )> v ( c fc> & fc)> v ( a fc,dfc)]- 

fe=0 

Proof. For proving equality (j2H|, we denote as /j and its left- and right-hand 
sides, respectively. We define also 

(7 fc = [z, k ,z k +v(a kl c k ) + v(&fc,dfc),v(cfc,dfc),v(afc,&fc)], 
ftfc = [zfc + v(a fc ,c fc ),z fc + v(&fc,dfc),v(cfc,6A!),v(ofc,dfc)], 

so that = nLo 3fc - IlLo h k- 

For n = 0, equality (|2"3)l is reduced to the identity Q22JI. Suppose that is 
true for n = 0, . . . , N — 1, N > 1. Then we have by induction 

= h N f^- 1] 

N—l 

+ [z N + v(b N ,c N ),z N + v(a N ,d N ),v(a N ,c N ),v(b N ,d N )] JJ 

JV-l JV 

= ?\ Yi '> k n ^ 

fc=0 fe=0 

where 

€n = h N + [zjv + v(6 w ,CAr),zjv + v(aAr ! dAr),v(ajv,CAr),v(6jv,djv)]- 
Using the Fay identity we find that ^ w = 9n- Therefore, = JlfcloSfc ~ 
IlfeLo ^ fc = /r j that is formula l|23|l is valid for arbitrary n. □ 

Remark 1. As noted by the referee, this theorem is a special case of the general 
construction for telescoping sums: 

n k—Xn n n 

fc=0 j=0 j=k+l 3=0 j=0 

whose proof follows similar lines: 

n k—1 n n k n n k— In 

j2( xk - f*) n x i n ^=5^11^ n 

fc=o i=o j=fe+i fe=oj=0 j=fc+i fe=oj=o j=fc 

with an obvious cancellation of extra terms in the latter sums. 
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For elliptic curves (i.e., for g — 1) formula (|23(l was proven by Warnaar |18j . 
Its further degeneration to the trigonometric level leads to a Macdonald identity, 
which was published for the first time in the Bhatnagar-Milne paper ,1 and which 
generalizes relations obtained by Chu in [2]. As shown in equalities of such 

type with special choices of parameters can be used for derivations of more fine 
structured summation and transformation formulas for bibasic and elliptic hyper- 
geometric series. 

We suppose that for some N > the points ajvj kjv> cjv> djv S S are such that we 
hit a zero of a theta function: [v(ajv, 6jv)] = or [v(cjv, cJjv)] = 0. Then, using the 
antisymmetry [v(cfc,6fc)] = — [v(6fc,Cfc)] equality (|23|l can be rewritten as 

[zfc + v(6 fc , c fc ), z fc + v(q fc , d fc ), v(afc, c fc ), v(fr fc , d fc )] 
^ [z fe + v(a fc ,c fe ),z fc + v(6fc,d/ c ),v(& fc ,c fe ),v(a fc ,d fe )] 

r Zj, Zj + vjcij, Cj) + v(bj, dj), \(c J ,d J ),\(a j ,b J )] 



E 



X 

Equivalently, for [v(co, bo)] = or [v(ao, do)] — we obtain the sum 
zfe + v(6fc, c^.), zfc + v(ofc, 4), v(ofc, Cfc), v(6fc, d fe )] 



,_ o [z ,z + v(a ,c ) + v(6 ,rfo),v(c ,(io),v(ao,6o)] 

x TT [zj, Zj- + vjaj, Cj) + v(frj, dj), v(cj, dj), v(aj, bj)] 

M I z i+i + v (a3+i,Cj+i),Zj+i +v(6j + i,dj + i),v(cj + i,6j + i),v(aj + i,dj+i)] 

_ TT [zfc,z fc +v(afc,Cfc) + v(&fc,dfc),v(cfc,d fc ),v(afc,frfc)] 
fcVi [ z fc + v (afe,Cfc),Zfc + v(fefc,dfc),v(cfc,fefc),v(a fe ,dfc)] ' 

In the elliptic case and its further degenerations, such equalities are useful in 
searching matrices whose inversions are given by simple analytical expressions. It 
is natural to expect that our relations will lead to a generalization of some of 
the Warnaar's results on multibasic hypergeometric sums and matrix inversions 
|18j . Moreover, it is worth to analyze possible extensions of the Krattenthaler's 
determinant from the Warnaar's elliptic case j 1 81 to arbitrary Riemann surfaces. 

By special choices of parameters we can try to give to the derived sums hyper- 
geometric type forms. We consider only two particular cases. Substituting 

z fc = u + v(xk), b k = c = x , Cfc ee x k , dfc = d , 



where u G C 9 , v(acfc) = v(xo,Xf.), k = 0,1,..., and u 2 = v(do,Xo), u 
v(xfc,afc) into (|26|l . we obtain the following identity. 

Corollary 2. 

^ [u + 2v(aifc)] [uq - u[ k) - u 2 , u[ k) ] ( [u + v(xj)] 



(fc) 



fc=0 



[uo] [uo-uf-ua.u^l+iV [v( ^ +l)] 

[u - u[ j) - U 2 + v(tfj),U^ + v(£j),U 2 +v(Xj)] 



[u^ +1) + U 2 + v(Xj+l),Uo - u^ +1) + v(xj+i),u - U 2 + v(Xj + i)] 
tt [up + v(xfc), up - u[ k) - u 2 + v(xfc), + v(a:fc), u 2 + v(x fc )] 
fc=i [v(acfc), u[ k) + u 2 + v(xfc), u - u[ k) + v(x k ),u Q - u 2 + v(x fe )] 



s 
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If there would exist a sequence of points a k such that u\ = , then four 
theta functions of the second factor in the left-hand side of (|27[1 would cancel 
each other and we would obtain an exact g > 1 analog of the sum l|15fl . Such a 
condition was assumed implicitly in the derivation of the corresponding formula 
in |15|. However, in general this is possible only for g = 1 (with the choice uj = 
du, a k = k + a,Q, x k = k + xq). As a result, our g > 1 analog of l|15Jl (which corrects 
the related formula of jTH]) does not obey all structural properties of the g — 1 
summation formula. Still, we draw attention to the fact that the right-hand side of 
(|27|l satisfies analogs of the balancing and well-poisedness conditions: 1) the sums 
of arguments of theta functions in the numerator and denominator are equal to 
each other; 2) the reciprocal theta functions in the numerator and denominator 
have arguments whose sums are equal to Uo + 2v(x k ). 

Another curious summation formula (inspired by a referee's suggestion) is ob- 
tained after substituting a k — a, c k = bo = c, b k+ \ = d k = x k for k = 0, 1, . . . , and 
Zk+i = z + v(c, x k ) for k > into 

Corollary 3. 



It is easy also to see that all the derived sums represent "totally abelian" func- 
tions, that is they are invariant under arbitrary moves of points on the Riemann 
surface along the cycles and appropriate 2g shifts of the variables z k (or uo). Anal- 
ogously, the Sp(2g, Z) modular group invariance is evident by construction (due 
to the Fay identity properties). As it is clear from our preliminary analysis, theta 
hypergeometric series for Riemann surfaces with the genus g > 1 should obey some 
principally new features with respect to the elliptic case and their determination 
requires some additional effort. 

The author is deeply indebted to the organizers of the Workshop on Jack, Hall- 
Littlewood and Macdonald polynomials (Edinburgh, September 23-26, 2003) for 
invitation to speak there and to the referee for pointing to a flaw in the original 
formulation of Corollary 2. This work is supported in part by the Russian Foun- 
dation for Basic Research (grant No. 03-01-00780). Its final version was prepared 
during the author's stay at the Max-Planck-Institut fur Mathematik in Bonn whose 
hospitality is gratefully acknowledged. 
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n 



[z ,z + v(c,x k -i) + v(a, x k ),-v(a,c),v(x k - 1 ,x k )] 




[z + v(a,x n ),v(c,x n )] _ [z + v(a,x ), v(c,x )] 
[z + v(c, x n ), v(a, x n )} [z + v(c, x ), v(a, x )] 
[z , z + v(c, xq) + v(a, x n ), v(a, c),v(x ,x n )] 
[z + v(c, x ), z + v(c, x n ), v(a, xq), v(a, x n )] ' 



(28) 
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